Levy-flight Spreading of Epidemic Processes leading to Percolating Clusters 
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We consider two stochastic processes, the Gribov process and the general epidemic process, that 
describe the spreading of an infectious disease. In contrast to the usually assumed case of short- 
range infections that lead, at the critical point, to directed and isotropic percolation respectively, 
we consider long-range infections with a probability distribution decaying in d dimensions with 
the distance as l/R d+rT . By means of Wilson's momentum shell renormalization-group recursion 
relations, the critical exponents characterizing the growing fractal clusters are calculated to first 
order in an e-expansion. It is shown that the long-range critical behavior changes continuously to 
its short-range counterpart for a decay exponent of the infection a — a c > 2. 
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I. INTRODUCTION 



A. Epidemic Processes 



The investigation of the formation and of the prop- 
erties of random structures has been an exciting topic 
in statistical physics for many years. In the case that 
the formation of such structures obeys local rules, these 
processes can often be expressed in the language of pop- 
ulation growth. It is well known that two special growth 
processes called (in the language of a disease) simple epi- 
demic with recovery (Gribov process also known in 
elementary particle physics as Reggeon field theory |g-§), 
the stochastic version of Schlogls first reaction and 
epidemic with removal (general epidemic process (GEP) 
P-|l0|) lead to random structures with the properties of 
percolation clusters: directed percolation |11-13| in the 
first case and isotropic percolation (for a recent overview 
see [0) in the last one |p^|-[l7|]. These stochastic pro- 
cesses describe the essential features of a vast number of 
growth phenomena of populations of infected individuals 
near their extinction threshold and are relevant to a wide 
range of models in physics, chemistry, biology, and soci- 
ology. The transition between survival and extinction of 
such a growing population is a nonequilibrium continu- 
ous phase transition phenomenon and is characterized by 
universal scaling laws. 

The Gribov process with short-range infection belongs 
to the universality class of local growth processes with 
absorbing states [p|Jl8| such as the contact process 1 19 2l}| 
and certain cellular automata |2^,^]. This universality 
class is characterized by the following four principles: 

1. Infection of susceptible ("birth") and spontaneous 
annihilation ("recovery") of infected individuals. 

2. Interaction ("saturation") between the infected in- 
dividuals. 



3. Diffusion ("spreading") of the disease in a d- 
dimensional space. 

4. The state without infected individuals is absorbing. 

To model these principles in a universal form, we use 
from the beginning a mesoscopic picture in which all mi- 
croscopic length- and time-scales are considered as very 
short. Thus we take a continuum approach with the den- 
sity n (x, t) of the infected individuals (the ills) at time 
t as an order parameter. Note that the spontaneous an- 
nihilation of the ills makes it possible to avoid compli- 
cations arising in the case of only diffusion controlled 
reactions which need creation and destruction operators 
as order parameters for a correct description. 

The Langevin equation is constructed in accordance 
with the four principles as 



d t n = W 2 n + R[n]n + (, 



(1) 



where £ (x, t) denotes a Gaussian Markovian noise with 
short range correlations which has to vanish if n (x, t) = 
to model the absorbing state and the reaction rate R [n] 
models birth, recovery and saturation. In a low density 
expansion we may set Va 



R[n] 



-At 



9 



(2) 



(C (x, t) C, (x', t')) = Xg'n (x, t) 6 (x - x') S (t - t') . (3) 

In contrast to the Gribov process (GP), the general 
epidemic process (GEP) introduces besides the suscep- 
tibles, S, who can catch the disease, and the infectives 
or ills, /, who have the disease and can transmit it, as a 
third class the removed, R, namely those who have had 
the disease and are now immune or death. Thus the first 
principle above is to be modified to 

1.' Infection of susceptible and spontaneous annihila- 
tion but without recovery of susceptible individu- 
als. 



1 



Therefore the reaction rate now also depends on the 
number of the removed individuals, introducing a mem- 
ory term into the process. Because this term is the lead- 
ing one in the long-time and large-distance limit we now 
have [d 



R [n] = —A (r + gm) , 
m(x,i) = A / dt'n(x,t') . 



(4) 
(5) 



In a microscopic realization of the Gribov process a sin- 
gle species of (quasi)particles, the Ps, is introduced. The 
J's represent the infected individuals (sites of a lattice). 
They perform simple random walks and undergo the fol- 
lowing "chemical" reaction scheme built from reversible 
branching and irreversible spontaneous annihilation: 



Gribov Process 



21 



Above some value of the branching rate of the reaction 
I — > 21, the stationary state has a finite density of J's. 
As the branching rate goes down to a threshold value, the 
stationary state density of J's goes continuously to zero, 
which is an absorbing state below the threshold. This 
threshold value corresponds to the critical point r = t c 
in Eq. ( |||). r c = as long as one neglects fluctuation 
contributions. 

A microscopic model which belongs to the universality 
class of the general epidemic process involves the three 
species S, I, and R specified above. Only the infected 
individuals / are mobile. A susceptible S may be contam- 
inated, but the infected Ps may become spontaneously 
immune: 



General Epidemic Process 



I - 
y R 
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The history of this model goes back to 1927 when it was 
first introduced in the mathematical biology literature 
|jTof . Here of course the stationary state is /-free. In 
this model the key parameter is the initial density of S"s, 
denoted by p: depending on the value p with respect to 
a threshold value p Cl the infected individuals may either 
start to proliferate as a solitary wave before dying out in a 
finite system, which occurs for p > p c (r < r c in Eq. (^)), 
or their number decreases from the outset, which takes 
place for p < p c (r > r c ). 



B. Levy-flight Infections 

In the standard version of the epidemic models the 
susceptible individuals can become contaminated by al- 
ready infected neighboring individuals. At the same time 
infected individuals are subject to spontaneous healing 
or immunization processes. 

In realistic situations the infection can be also long- 
ranged. This may be e.g. by a disease in an orchard where 



flying parasites contaminate the trees practically instan- 
taneous in a widespread manner if the timescale of the 
flights of the parasites is much shorter as the mesoscopic 
timescale of the epidemic process itself. Thus following a 
suggestion of Mollison || , Grassberger uM introduced a 
variation of the epidemic processes with infection proba- 
bility distributions P (R) which decays with the distance 
R as a power-law like 



P(R) oc 



1 



R d+o 



for R — ► oo. 



(6) 



We will in general denote such long-range distributions as 
Levy-flights although a true Levy-flight is defined by its 
Fourier transform as P (q) oc exp (—6 |q| ff ) 25|, and then 
only Levy-exponents with < a < 2 give rise to positive 
distributions [ p6| . The infection rate in the Langevin 
equation ([!]) is now given by 



dn (x, t) 



dt 



After Fourier transformation of this equation and after a 
small momentum expansion that is relevant in our meso- 
scopic consideration we get 



dh (q, t) 



dt 



(Po 



P2Q 2 + Vo<f + {q 2 , q a )) h (q, t) 



inf 



(8) 



where the analytical terms stem from the short-range 
part of P (R) and the nonanalytical ones arise from the 
power-law decay (jfy. The constant po is included in the 
reaction rate as a negative ("birth") contribution to t 
whereas pi yields a diffusional term. Naively the param- 
eter pu is relevant or irrelevant in the long wave-length 
limit if a is smaller or bigger than 2 respectively, and this 
fact has mislead some authors to neglect this term from 
the outset if a > 2. But this naive ("Gaussian") argu- 
mentation may be wrong in an interacting theory because 
the critical behavior is in general determined by an non- 
trivial fixed point of a renormalization group transfor- 
mation. To decide which one of the terms in Eq. (||) are 
relevant, one has to compare with the scaling behavior 
of the Fourier transformed susceptibility \ (Qj w ) °^ q 2 ~ Tl ■ 
If a < 2 — fj, the parameter p a is a relevant perturbation 
and must be included in a renormalization group pro- 
cedure. Prominent examples of systems with fj < are 
</> 3 -models as e.g. the Yang-Lee-singularity model |27]]. In 
all these cases p a is relevant also for a > 2 and cannot 
be neglected. 

In the following we define a — 2 (1 — a) and the dif- 
fusion term in Eq. is now completed by a term pro- 
portional to ex g 2 ( 1-a )ro (q,t). In real space we write the 
completed Langevin equation as 



dm — A 



v 

2c\ 



-v 2 ) 



V 2 n + R [n] n + £ 



(9) 
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and the gradient-terms should be only considered (in 
Fourier space) up to a cutoff A that we have set to 1. 
Then stability of this terms against inhomogeneous per- 
turbation is guaranteed if v > 0. 

Grassberger 24 reported new critical exponents for 
a > from a 1-loop calculation that contain some nu- 
merical errors. These exponents are discontinuous in the 
limit a — > +0 if one assumes irrelevance of the new terms 
for a < 0. In this paper we will reconsider the problem 
and show that the full range of values a > fjsn/2 lead 
to new critical behavior. Here fjgR < is the anoma- 
lous susceptibility exponent of the epidemic models with 
short-range infection. We will show that the critical ex- 
ponents change continuously at the boundary 2a ~ fjsB. 
from long-range to short-range behavior. 

We remark that the interest in reaction-diffusion prob- 
lems involving particles that perform Levy flights is not 
new. In the physics literature, they have most recently 
arisen as follows. Particles performing simple random 
walks are subject to the reactions A + B — > and 
A + A — > in the presence of a quenched velocity field 
p8j. The effect of the quenched velocity field is then 
to enhance diffusion in such a way that the effective ac- 
tion of the velocity field is reproduced if Levy flights are 
substituted for the simple random walk motion. In the 
above mentioned reactions the time decay of the particle 
density is algebraic with an exponent related to that of 
the step length distribution of the Levy flights defined 
in Eq. (|^) . These results have been confirmed by several 
renormalization group calculations [§9|,|o) . 



II. THE GRIBOV PROCESS WITH 
LEVY-FLIGHTS 

A. Renormalization Group Analysis 

In order to develop the renormalization group analysis 
we recast the Langevin equation (||) as a dynamic func- 
tional lEnT-" 



J[s,s] = J d d xdt~s\d t + \^{s-s) 



T- V 2 



+ V 2 



}s (10) 



where s is Martin-Siggia-Rose response field |Q. We 
note that he dynamic functional can also be derived us- 
ing the methods developed in [ 35|-p7t from a microscopic 
master equation. By a suitable rescaling of the density 
n oc s, the constant g' i n Eq. (||) is made equal to g. The 
dynamic functional (|l0|) is then symmetric in the absorb- 
ing phase under the exchange s {x,t) — s (x, — t). All 
correlation and response functions can be calculated as 
functional integrals with weight exp {—J) in a pertur- 
bation expansion involving the propagator (the unper- 
turbed response function) 



Go (q,i) = e(t)exp{- 



2a 



} 
(11) 



as a function of momentum q and time t. This prop- 
agator guarantees stability for all a as long as r > 0, 
v > 0, and q = |q| < 1. For simplicity we have set the 
momentum cut-off A = 1. 

To study the critical behavior of this system near the 
critical point we use Wilson's renormalization procedure. 
We introduce the usual coarse graining parameter b > 1 
and split the fields s and s into components which are non 
zero on the momentum shell fib = {q| 1 /fo < |q| < 1} and 
components defined on the complement of Qb, the latter 
being denoted by s K and s < . We integrate out the short 
scale degrees of freedom in the weight exp {—J), that is, 
those defined on f^, and rescale the fields according to 

s (x, t) - s' (b- 1 *, b- 2 -^t) - b^/ 2 s < (x, t) , 

S (X, t) -► S 1 (b- 1 *, b~ 2 ^t) = fo(d+T)/2 5 < ( Xj f J _ ( 12 ) 

Rcnormalizcd parameters r', v' , and g' arc defined in 
such a way that the coarse grained functional looks like 
the old one. The one-loop calculation is standard and 
does not present any technical difficulties. For infinites- 
imal renormalization transformation with b — 1 <C 1 wc 
obtain 



iuj + A 



;6" 7 



r' + q 2 + ^(q- 2a -l) q 2 



4(l+ry 



lnfe 



-\b 2+ ^ 



2 (l + r) 



In b 



+\q 2 b<-< 1 + i- (b 2a q- 2a - I) 



uK (v) 
8(1 + t) : 



In 6 



and 



2u 



(1+r 



lnfe 



(13) 



(14) 



Here u> is the frequency, u = Sag 2 /2 with Sd the surface 
of the unit sphere in d dimensions divided by (2w) , and 
K (v) = 1 — cv, where c is an uninteresting positive con- 
stant. Note that the calculation neglects terms of order 



O 



but is exact (for the coarse graining method we 



have used) with respect to the parameter v. By com- 
parison of the terms oc u and q 2 in Eq. (Jl|) we get the 
Wilson- functions 



7 = 



4 (i+7-r 



O 



7 = 7 - C 



uK (v) 
8 (l + r) 5 



O u 2 



(15) 
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We use I = In b as the flow parameter of the renormal- 
ization transformation that yields then from the other 
terms in Eqs. dl3,n4) the flow equations 



dr 
~dl 

du 
~dl 



(2-7)r + 



2(l + r) 



+ O (u 2 ) 



4 - d — 7 - 27 ■ 



2 a 



- = (2a- 7 ) U . 



(16) 
(17) 
(18) 



The last equation is exact sinct 
-V 2 ) a s in the dynamic functional 



the operator 

s [—v j sm tne dynamic tunctionai ^ |i0| ) is not renor- 
malized. The reason is that the renormalization proce- 
dure can only generate contributions that are analytic in 
the momenta. 

In order to study the fixed point structure we find it 
useful to introduce £ = 4 (1 — a) — d = e — 4a. Near the 
fixed points we linearize the flow equation ( |l6| ) for the rel- 
evant variable r about t c = r* as dr/dl ~ v~ x (r — r*). 
The flow equations ( p|Jilj| ) have, besides the trivial short- 
range Gaussian fixed point (t*,u*,v*) = (0,0,0), stable 
for a < and d > 4, and the trivial Levy-Gaussian 
fixed point (r*,u*,v*) = (0,0,2a) with t)lg = 7* = 2a, 
z = 2 + C* = 2 (1 - a), stable for d > 4(1 - a) and 
a > 0, two non trivial fixed points. The first one is the al- 
ready known short-range directed percolation fixed point 
(t*,u*,v*) = (-£/3,2e/3,0) + O(e 2 ), with fj DP = 7* = 
-e/12+O (e 2 ) , zup = 2 + C* = 2-e/12 + (e 2 ) , v DP = 
2-e/4 + (e 2 ) , stable for e > 0, a < -e/24 + O (e 2 ) 
. The second one is the new Levy-directed-percolation 
fixed point 



28a + e 



14- 



O {e 2 ) 

£ 



7(1 - a) 



+ O e 2 



(19) 



which is stable for -e/24 + (e 2 ) < a < e/4. We obtain 
for this fixed point the critical exponents 



Vldp = 7* = -y + 0(e 2 ) , 
?7ldp = 7* = 2a, 

Z LDP = 2 + ?|iDP - fjLDP 

= 2{l-a)- £ -+0 (e 2 ) , 
^i,P = 2(l-«)-f + 0(e 2 ). 



(20) 
(21) 

(22) 
(23) 



Note that at all fixed points u* and v* are non-negative 
as they should for stability of the theory. 

We have depicted in Fig. |l] the stability regions for 
each of the above fixed points in the (a, d) plane. A 
glance on the exact flow equation (|lj) of the parameter 



v shows that the boundary between the domains of at- 
traction of the directed percolation fixed point and the 
Levy-directed-percolation fixed point is given exactly by 
fjop = fjLDP = 2a. At this boundary all exponents 
change continuously upon varying the parameter a from 
their values at the directed percolation fixed point to 
those of the Levy-directed-percolation fixed point and 
vice versa. 

d-4 



DP 



-24a 



LG 



> a 



4a 



LDP 



FIG. 1. Stability regions of the Gribov process with long 
range spreading in the (a, d — 4 — e) plane. G, LG, DP 
and LDP indicate the respective stability regions of the short 
range Gaussian, Levy Gaussian, short range directed perco- 
lation and Levy directed percolation fixed point 



B. Scaling analysis 

In the following we consider the scaling behavior of 
two key quantities: the time dependent order parameter 
(the density of infected individuals) p(t) — (s(x, i)) Po 
for t > if the initial state at time t = is prepared 
with a homogeneous initial density po , and the response 
function x(x, t) = (s (x, t) s (0, 0)) that yields the density 
of infected individuals after the epidemic is initialized 
by a pointlike source at t = and x = 0. Here we 
are interested in the Levy-flight case only, thus we will 
disregard the subscripts at all critical exponents. We set 
t — t c — ► t in the following. 

^From the rescaling Eqs. (|12|) we get at the Levy fixed 
point (which is approached for b 3> 1) the relationship 

p (t, r, p ) = b-^' 2 p (b-% bV"T, b^ 2 P0 ) , (24) 

where the critical exponents are displayed by Eqs. ( p0| - 
^3| ) . The scaling of the initial density is easily found by 
noting that po arises in the Langevin equation (^J) as an 
additive source term g(x, t) = pg5 (t) that translates in 
the dynamic functional J, Eq. (]10|), to a further addi- 
tive contribution J d d x pos(x, 0) from which one directly 
reads off the scaling behavior of po if one knows that 
s(x, 0) scales as s(x, t) |38[], At criticality, for r = 0, 



4 



one obtains from Eq. (B4J) that the order parameter first 
increases in a universal initial time regime J3^,|33| as 

p(t,p )xp t e (25) 

where the universal initial time exponent 9 is given by 

a _ V _ e 



Z 7(7 



+o(. 2 ), 



(26) 



and we have set a = 2(1 — a). Then, after some 
crossover time, the order parameter decreases as 



with 



d + rj „ 3e 



2z 



(27) 



In the non absorbing stationary state, i.e. for t < and 
t — > oo, the order parameter behaves as 

Pstat (r) cx Irl' 3 , (28) 
where the order parameter exponent (3 is found as 
^(d + 77) 



(29) 



Finally, at the critical dimension rf = 2cr = 4(1 — a), 
the scaling behavior is mean-field like with logarithmic 
corrections, 

p(t) cx p \n 1/7 t 

in the initial time regime at criticality, 

PW K — — 

in the long time regime at criticality, 

p stat (r) cx|r|ln 2 / 7 (^\ 

in the stationary state, (30) 

which we mention for completeness. 

The scaling behavior of the response function is given 

by 



X (x, t) = b~ {d+7>) x (V'x, b~% fcVO 



(31) 



First we read off the correlation lengths for space and 
time and the corresponding exponents as 



with 

£|| « l T 
with i/ii 

At criticality, r = 0, we have 



v± = v = - + t^+O (e 2 



zv = 1 + — +0 (e 2 ) 



(32) 
(33) 

(34) 



with a universal scaling function T (x) . This relation 
shows that the density of infected individuals as the re- 
sult from a pointlike seed dies out with an exponent 
(d + rj)/ z — 2(3/ uz. Comparing with the general decay 
law d27n , we find that the probability to find after a time 
t an infected individual if there was a pointlike seed at 
time t = decays with an exponent (3/vz. The Fourier 
transformed susceptibility at the critical point scales as 



X (ci,u;) = q- 2+11 T(cj/q z ) X q- 



(35) 



As the last scaling exponent that can be deduced in the 
usual way from the given critical exponents we present 
the fractal dimension of the clusters of the infected indi- 
viduals: 



df 



V 



d-? = 

v 2 



(36) 



We note that the values of the exponents v\\ (J32J), v\_ (|33j), 
j3 (^9|) , and d / (|3^) that we have found are different from 
those given by Grassberger |2^| . The values of all the 
exponents changes continuously at the stability boundary 
cr = 2 (1- a) = 2- fj DP = 2 + e/12 + O (e 2 ) to their 
short-range directed percolation values. 



C. Comparison with existing simulations 

In a recent letter Albano Q] has presented a numer- 
ical study of one-dimensional branching and annihilat- 
ing random walks (BARW) in which the individuals per- 
form Levy flights. For Brownian particles in d < 2 di- 
mensions the BARW (which is defined by the equations 
A — > (to + I) A and A + A — > 0) is known to belong 
for to odd to the universality class of directed percola- 
tion . For d > 2 the systems shows a phase tran- 
sition at zero branching rate which can be described by 
mean field exponents. If the random walk is replaced by 
Levy flights, noise becomes irrelevant above d a = a. 

Albano has investigated the behavior of the critical 
exponents for m = 1 as a function of ct, for 0.25 < a < 11. 
His results are summarized in the following table: 



a 


z 


V 


z — a — rj 


2 


1.590 


-0.482 


0.072 


1.5 


1.585 


-0.483 


0.568 


1 


1.583 


-0.489 


1.073 


0.75 


1.581 


-0.512 


1.343 


0.5 


1.575 


-0.553 


1.628 


0.25 


1.569 


-0.574 


1.893 



(37) 



The critical dimension d a is lower than 1 for a G 
{0.25,0.5,0.75}. (For a = 0.25 even the critical di- 
mension of the Levy-flight directed percolation d c — 
4(1 — a) = 2<7 is lower than 1). Therefore the phase 



5 



transition should occur at zero branching rate, and crit- 
ical exponents should be the mean field ones (/3 m f = 1, 
z m t = & and r) m f = 0); this is clearly not the case for the 
exponents given in table (|37f). The hyperscaling relation 
z = a + r/ Eqs. ( |2l| , |22| ), which we have shown to hold 
in any dimension, is violated as a is decreased. These 
facts cast doubt on the reliability of the simulations per- 
formed in |40(| . A possible explanation is to be found in 
the Levy flight generation procedure. Indeed the author 
uses a distribution Eq. (||) truncated at some distance 
cut-off, the effect of which is to produce an effectively 
short range motion. This interpretation is confirmed by 
the slow variation of the exponents as a function of <r, 
their values remaining close to that of directed percola- 
tion with simple random walk displacement. 



III. THE GENERAL EPIDEMIC PROCESS WITH 
LEVY-FLIGHTS 

A. Renormalization Group Analysis 

The renormalization group analysis of the GEP is per- 
formed analogously to the corresponding analysis of the 
Gribov process presented in section 2. The Langevin 
equation for the GEP @ , where now the reaction rate is 
given by Eq. (JI]) , is recast in the dynamic functional [jl6| : 

J [3,a] = J d d xdt~s\d t + XgS - A|s 



+A 



)]}'■ (38) 



The field S (x, t) — A dt' s (x, t'), a rescaled form of 
Eq. (|| ), introduces a memory term in the dynamics. In 
analogy to the Gribov process we have rescaled the fields 
so that g 1 = g. The dynamic functional (38) is then 
symmetric under the exchange S (x, — t) <-> — s (x, t), or 
s (x, — t) «-> <9ts(x, t) |16[] . From this symmetry follows 
that we only need to consider one coupling coefficient 
g for the two interaction terms in J . The perturba- 
tion expansion involves also the propagator displayed in 
Eq. (p"l|), and v > is needed for stability. We inte- 
grate out the short scale degrees of freedom in the weight 
exp (-J), and rescale now the fields according to 

s (x, t) -> a' (6- x x, b- 2 -<t) = 6 ^+2+7)/2 s < (Xj t) t 

§ (x, t) - s' (6- x x, b- 2 -<t) = &( d - 2 +W 2 S< ( Xj t ) . (39) 

Note that from the exchange symmetry follows exactly 



c = 



7-7 



(40) 



The renormalized parameters r', v', and g' are now de- 
fined by the coarse graining equations calculated to one- 
loop order with b — 1 <C 1 



2a 



zcj5 -(7+7)/2 
+ \b 2 -~~< 

+X q 2 b~^ 



I - 



3u 



4(1 + t) 
In b 



■Info 



2(1 + tY 



la -la 



— (b 2 °q 
2a 1 y 

(d - 2) uK (v) 

4d(l + r) 2 



1 



In b 



and 



ub 6 ~ d -^ 



I - 



4?/ 



(1 



■In 6 



(41) 



(42) 



where u = Sdg 2 , and K (v, r = 0) = 1 — cu, and c is an 
uninteresting positive constant. The Wilson functions 
and the renormalization group equations follows now as 



7 + 7 



dr 

~dT 

du 

~dl 

dv 

~dl 



3u 



2(1+T) 

(d - 2) uK (v 
4d(l + 

(2 - 7) t + 



O (u 2 ) , 
— + (« 2 ) . 



a 



2(1 + T y 

4h 



O u z 



6 - d - 3 7 ^— r + O (u 2 ) 

(1 + r) 3 V ; 

(2a-7)u, 



u, 



(43) 
(44) 

(45) 

(46) 
(47) 



The last equation is exact. 

In order to study the fixed point structure we find it 
useful here to introduce e = 6(1 — a) — d = e — 6a. 
The flow equations ( f45f]47| ) have, besides the trivial short- 
range Gaussian fixed point (t*,u*,v*) — (0,0,0), stable 
for a < and d > 6, and the trivial Levy-Gaussian fixed 
point (t*,u*,v*) = (0,0,2a) with t)lg = 7* = 2a, t] LG = 
7* = -2a, z = 2 + C* = 2 (1 - a), stable for d > 6 (1 - a) 
and a > 0, two non trivial fixed points. The first one is 
the already known short-range- GEP fixed point (where 
the static exponents are known from undirected percola- 
tion) (t*,u*,v*) = (-e/7,2e/7,0)+O (e 2 ), with fj GEP = 
7* = s/21 + O (e 2 ) , r lGEP = 7 * = -8e/21 + O (e 2 ) , 
z GEP = 2 + C* = 2 - e/6 + O (e 2 ) , ^ep = 2 - 5e/21 + 
O (e 2 ), stable for e > 0, a < -e/42 + O (e 2 ). The second 
one is the new Levy-GEP fixed point 

u* = le + 0(e 2 ), 



96a(l-a) + (2-3a)f? 



48(1 - a) +c (2 

e 



T =T r = - 



16(1 - a) 



3a) e 
h O (e 2 



+ o ( e - 2 ) 



(48) 
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which is stable for — e/42 + O (e 2 ) < a < e/6. We obtain 
for this fixed point the critical exponents 



Vlgep = 7* = -2" - tt + ( £2 ) > 
o 

Vlgep — 1* = 2a, 

z lgep = 2 + C = 2 (1 - a) - j| + O (e 2 ) , 



"Zgsp = 2(l-a)-J + (e 2 ) 



(49) 
(50) 
(51) 

(52) 



Note again that at all fixed points w* and v* are non- 
negative as they should for stability of the theory. We 
have depicted in Fig.|| the stability regions for each of 
the above fixed points in the (a, d) plane. Now a glance 
at the exact flow equation ( |47| ) of the parameter v shows 
that the boundary between the domains of attraction of 
the short-range-GEP fixed point and the Levy-GEP fixed 
point is given exactly by fjcEP — Vlgep — 2a. At this 
boundary all exponents are again continuous functions of 
the parameter a. 



d-e 



G 



GEP 



e = -42a ! 



LG 



-> a 



e = 6a 



LGEP 



FIG. 2. Stability regions of the general epidemic process 
with long range spreading in the (a, d = 6 — e) plane. G, LG, 
GEP and LGEP represent the stability regions of the short 
range Gaussian, Levy Gaussian, short range percolation and 
Levy percolation fixed point, respectively. 



B. Scaling analysis 

At first we consider the scaling behavior of the time de- 
pendent density of infected individuals p (t) = (s (x, t)) Po 
for t > if the initial state at time t = is prepared with 
a homogeneous initial density po- We are interested in 
the Levy-flight case only, thus we will disregard the sub- 
scripts at all critical exponents. Again we set r — r c — > r. 

Here we find from the rescaling Eqs. ( |39| ) at the Levy 
fixed point (which is approached for b » 1) the relation- 
ship 

p (t, T, po) - h -id+^)/2 p U- % b l/v T} b (d+2-n)/2 pQ J ^ 



(53) 



where the critical exponents are displayed by Eqs. ( [49| - 
|52|) . The scaling of the initial density is also found by 
adding a source term (/(x, t) — po5 (t) to the Langevin 
equation (p|) that translates in the dynamic functional 
J , Eq. (pq), to the additive contribution J d d x pos(x, 0) 
from which one reads off the scaling behavior of po since 
s(x, 0) scales as s(x, t). At criticality, for r = 0, it follows 
from Eq. ( |53| ) that the infection density first increases in 
the universal initial time regime as 

p(t,po) tx p t e (54) 

where the universal initial time exponent 6 is given by 

V + V 



2z 



o~ 

z 



1 



3s 
Tfkr 



+0 ( £ - 2 ) 



(55) 



We have set a = 2 (1 — a). 

As an order parameter we consider the density of the 
removed (immune) individuals (^|) namely 



p(t,T,p ) oc (S (x,t)) f 



= A 



df (s (x,i)) P0 . 



(56) 



The scaling properties of this order parameter are deter- 
mined from Eq. (|5^) by 

p (t, r, po) - 6-( d - 2 +^)/ 2 p (b- z t, b^r, b^ d+2 -^/ 2 po) . 

(57) 

The initial infection density po pl a Y s here the role of an 
ordering field. In the infinite time limit at criticality, 
when t = 0, the order parameter p sta t = p(t — > oo) goes 
to zero with po as 



Pstat (t = 0,po) oc p Q 1/S , 

~ d + 2-fj 



with 



d-2 



V 



(58) 



Below threshold, that is r > 0, the order parameter is 
linear in po with a coefficient that diverges as r — > + : 



Pstat (t > 0,p ) oc p T 7 , 
with 7 = (2 — fj) v 

e 



1 



4(7 



+0 e 



2^ 



(59) 



Lastly for r < 0, the order parameter is independent of 
Po in the limit po —> 0, and goes to zero with r as 

Pstat (t < 0) OC |^|^ , (60) 

where the order parameter exponent is given by 
(d-2- 



■fj) (d- 
— = v — 



1 - ^+0 (e 2 ) 
4cr v ' 



(61) 
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At the critical dimension d — 3a, the scaling behavior of 

Pstat 

pstat(T = 0, pa) oc p 1/2 

ln^r 



yO s tat(T > 0, pa) OC p - 



Pstat(T < 0) oc |r| In 



1/4 



(62) 



is mean-field like with logarithmic corrections. At criti- 
cality we find from (|57|) the scaling behavior 



(63) 



Pl^PoJ — Po f [tPo 

with a universal scaling function 

t"l z for x < 1 

(exponentially) for x — > oo. 



J~ (x) oc -j 1 



(64) 



Note that the "static" exponents /3, 7, J, 77 (and j/) cor- 
respond to well known undirected percolation exponents 
but for long-range connectivity. They were already given 

in US 

To study the spread of the infection by computer 
simulations one may investigate the response function 
X(x, t) = {S (x, t) s(0, 0)) which describes the density of 
the immune percolating individuals at the time t caused 
by an infection at t = 0, x = 0. Its scaling behavior is 
given by 

X(x, t, t) = b-^+^x (6 _1 x, b- z t, 6 1/l/ r) . (65) 

We read off the correlation lengths for space and time 
and the corresponding exponents as 



1 £ 

with v = - + —^+0 (e 2 ) , 
a 4eH 

with v t = zv = 1 + — — hO (e 2 ) . 

16(7 

At criticality, r = 0, we have 

x(x,t) = |xr^ 2+f ')j- x (|x|/t^) 



(66) 
(67) 

(68) 



with an universal scaling function T (x) . The Fourier 
transformed susceptibility at the critical point scales as 



X (q, W ) = (? - 2 +''- 2 ^(^)oc 



(69) 



As the last scaling exponent that can be deduced in the 
usual way from the given critical exponents we present 
the fractal dimension of the percolation clusters of the 
removed individuals: 



df = d = — 

v 



(70) 



We note that the value of the exponent vt (|67|) that we 
have found is different from that given by Grassberger 
|p4[| . The values of all the exponents change continuously 
at the stability boundary cr = 2 (1 — a) = 2 — fjcEP = 2 + 
e/21 + O (e 2 ) to their short-range undirected percolation 
values. 



IV. CONCLUSIONS 

Epidemic processes are growth models for phenom- 
ena arising abundantly in nature. We have shown by 
imposing a Levy flight type of infection spreading that 
new long-range determined universality classes come into 
play. We were able to characterize the universality classes 
by determining the critical exponents to first order in 
an e-expansion around their upper critical dimension. 
There exist exact relationships between the exponents, 
and some critical exponents are given exactly as functions 
of spatial dimension and the exponent characterizing the 
long-range tail of the Levy flight infection. Besides, we 
have been able to build a renormalization group flow that 
possess a fixed point structure that allows to describe 
short-range and long-range infection in the same and to 
pass continuously from one behavior to another by vary- 
ing the Levy flight exponent. Because the anomalous 
susceptibility exponent, the analog of the Fisher expo- 
nent in critical equilibrium phenomena, is negative here, 
the continuous crossover between long- and short-range 
behavior arises at a Levy exponent greater than 2. We 
hope that this work triggers more simulational work on 
this subject. 
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